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ABSTRACT 

We describe the operation and performance of the difference imaging pipeline (Diffimg) used to 
detect transients in deep images from the Dark Energy Survey Supernova program (DES-SN) in its 
first observing season from 2013 August through 2014 February. DES-SN is a search for transients in 
which ten 3-deg^ fields are repeatedly observed in the g,r,i,z passbands with a cadence of about 1 
week. The observing strategy has been optimized to measure high-quality light curves and redshifts 
for thousands of Type la supernova (SN la) with the goal of measuring dark energy parameters. The 
essential Diffimg functions are to align each search image to a deep reference image, do a pixel- 
by-pixel subtraction, and then examine the subtracted image for significant positive detections of 
point-source objects. The vast majority of detections are subtraction artifacts, but after selection 
requirements and image filtering with an automated scanning program, there are ~ 130 detections 
per deg® per observation in each band, of which only ^ 25% are artifacts. Of the ^ 7500 transients 
discovered by DES-SN in its first observing season, each requiring a detection on at least two separate 
nights, Monte Carlo (MG) simulations predict that 27% are expected to be SNe la or core-collapse 
SNe. Another ~ 30% of the transients are artifacts in which a small number of observations satisfy the 
selection criteria for a single-epoch detection. Spectroscopic analysis shows that most of the remaining 
transients are AGN and variable stars. Fake SNe la are overlaid onto the images to rigorously evaluate 
detection efficiencies and to understand the Diffimg performance. The Diffimg efficiency measured 
with fake SNe agrees well with expectations from a MG simulation that uses analytical calculations 
of the fluxes and their uncertainties. In our 8 “shallow” fields with single-epoch 50% completeness 
depth ~ 23.5, the SN la efficiency falls to 1/2 at redshift z ss 0.7; in our 2 “deep” fields with mag- 
depth ~ 24.5, the efficiency falls to 1/2 at zRi 1.1. A remaining performance issue is that the 
measured fluxes have additional scatter (beyond Poisson fluctuations) that increases with the host 
galaxy surface brightness at the transient location. This bright-galaxy issue has minimal impact on 
the SNe la program,but it may lower the efficiency for finding fainter transients on bright galaxies. 

Subject headings: techniques: image processing, supernovae 
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1. INTRODUCTION 

The di scovery of the accelerating expansion of the 
universe (jRiess et al.l 119981 : iPerlmutter et alJ Il999[ i us¬ 
ing Type la supernovae (SN la) has greatly motivated 
ever larger transient searches in broadband imaging sur¬ 
veys. The associated search pipelines have become in¬ 
creasingly complex in distributing enormous computing 
tasks needed to rapidly find new transients for spectro¬ 
scopic observations, and in processing a wide range of 
data quality. 

A new era of transient searches began in the early 2000s 
with “rolling searches” in which the same telescope is 
used for discovering new objects and providing precise 
photometric measurements of the light curve in multiple 
passbands. To collect large SN la samples for measuring 
cosmological parameters, the earliest rolling searches in¬ 
clude the Supernova Lega cy Survey fSN LS: lAstier et ^ 


2006 : IPerrett et ^ 1201011 . ESSENCE (iMiknaitis et al 


200711. and the Sloan Digital Sky S urvey-II (SDSS-II: 
Frieman et al.ll200^ iSako et mi2008[ l . Each of these sur¬ 


veys discovered many hundreds of SNe la, about half 
of which were spectroscopically confirmed. The next 
generation of rolling se arches includes the recently com¬ 
pleted Pan-STARRSl (iKaiser et al.l 1200211. the on going 
Dark Energy Survey IDES: iBernstein et al. 12012117 and 
the Large Synoptic Survey Telescope (LSST: Ivezic et 'HI 


120081: IlSST Science Collaboration et al.ll2009l l . expected 
to begin in the next decade. Another advantage of these 
rolling searches is that there is a complementary wide- 
area survey with the same instrument; this benefits the 
absolute calibration by including dithered exposures over 
the SN fields to inter calibrate the CCDs, regular obser¬ 
vations of standard-star fields, and measurements of the 
telescope and atmospheric transmission functions. 

The goal of this paper is to describe the difference¬ 
imaging pipeline (Diff Img) used to discover point-source 
transients in DES. We present detailed performance re¬ 
sults of Diff Img for single-epoch detections, and for the 
redshift dependence of discovering and classifying SN la 
light curves. While the search strategy was optimized to 
find SNe la to build a Hubble diagram for measuring dark 
energy properties, (jBernstein et al.ll2012ll . Diff Img does 
not depend on the transient type. In addition to SNe la, 
our Diffimg has found many other transient types in¬ 
cluding c ore-collapse SNe (CC SNe) , super-luminous SNe 
(SLSNe: iPanadoDOulos et ^I2ni5h . activ e galactic nu¬ 
clei ( AGNs), Kuiper belt objects (KBOs: iGerdes et alJ 
201511. and a pos sible tidal disruption event ( Reeslll988t 


Folev et aI1l26l5ll . 


The challenges for Diffimg are to produce a high qual¬ 
ity subtracted image for each search image by subtract¬ 
ing a deep coadded template, reject a large number of 
non-astrophysical detections (artifacts) in the subtracted 
images, develop a workflow to process each night of data 
in less than a day, and monitor the performance well 
enough to uncover subtle problems and to determine ef¬ 
ficiencies and biases for science analyses. We use publicly 
available codes for the core routines needed to determine 
an astrometric solution, co-add exposures, measure the 
point-spread function (PSF), align template and search 
images, perform the subtractions, and fit light curves 
to a series of SN templates for classification. In addi¬ 
tion to these existing codes, we have developed new soft- 


ware tools for automa ted scanning of subtracted images 
(|Goldstein et al.l[2015L hereafter G15), detailed monitor¬ 
ing based on artificial SNe overlaid on images, and a 
workflow to distribute jobs on arbitrary computing plat¬ 
forms. 

The essential monitoring element is to inject fake 
SNe la onto galaxies in real images (hereafter called 
“fakes”). The Supernova Gosmology Project used fakes 
to monitor the effic iency of human s canners in the real¬ 
time SN la search (jPain et al.ll2002tl . and also to mea¬ 
sure the analysis efficiency as part of the SN la rate 
measurement. Fakes were later used for re al-time mon¬ 
itorin g in the SDSS-II Supernova search (|Dildav et al.l 
1200811 to measure the efficiency of the detection pipeline 
and human scanning. The Nearby Supernova Factory 
moved stars on the image to serve as fake transients; 
they monitored their single-epoch detection efficiency 
and trained their machine learning method that was 
used to reject lar ge numbers of subtraction artifacts 
(iBailev et al.ll2007H . SNLS used fakes in an offline anal¬ 


ysis ( Perrett et al.ll20in hereafter PIO) to measure their 
efficiencies and selection biases that impact the Hubble 
diagram. In DES-SN the fakes are used to (I) mon¬ 
itor the detection depth, (2) monitor the single-epoch 
detection efficiency from the fraction of fakes that are 
detected, (3) monitor the efficiency for multiple detec¬ 
tions that are required for spectroscopic targeting and 
science analysis, (4) train the automated scanning soft¬ 
ware (G15), and (5) characterize the Diffimg perfor¬ 
mance for a fast Monte Carlo (MG) simulation. 

Compared to the use of fakes in previous surveys, an 
improvement in DES-SN is that the ideal efficiency is 
predicted from the signal-to-noise ratio (S/N) of the flux 
measurements, and thus the Diffimg performance can 
be rigorously evaluated by comparing the predicted and 
fake efficiencies. This prediction, as a function of fake SN 
redshift, comes from a fast MC simulation that computes 
realistic light curves without using images or pixels. The 
fast MC simulation analytically computes the light curve 
fluxes and their uncertainties using input from observed 
conditions, and also from key Diffimg properties derived 
from the fakes. The agreement (or lack of) between the 
predicted and measured efficiency provides a robust mea¬ 
sure of the Diffimg performance. 

There is another practical motivation for using a fast 
MC simulation to validate the point-source Diffimg effi¬ 
ciency with fakes. Typical science analyses require large 
SN simulations that are repeated many times for devel¬ 
opment, evaluation of systematic uncertainties, and es¬ 
timates of contamination from CC SNe. Ideally, such 
simulations would be similar to the fakes in which calcu¬ 
lated light curves are overlaid on CCD images and pro¬ 
cessed with Diffimg. The CPU resources for so many 
image-based simulations, however, would be quite enor¬ 
mous. On the other han d, the fast MC simulation in 
SNANA (|Kessler et al.ll200^ can generate close to 10^ light 
curves per second on a single core, which is five orders 
of magnitude faster than the ideal image-based simula¬ 
tion. Our goal, therefore, is to use a single realization 
of fakes to characterize the Diffimg performance for the 
fast MC simulation; the fast MC can then be used to 
rapidly generate samples of point-source transients with 
the same efficiencies and uncertainties as an image-based 
simulation. Although only one transient type (SN la) is 
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used to generate fakes for image overlays and Diffimg 
processing, the resulting fast MC simulation can in prin¬ 
ciple be used for any SN type, and more generally for 
any point-source transient. 

The outline of this paper is as follows. An overview of 
DES and the transient search is given in ^ and Diffimg 
is described in The monitoring of single-epoch detec¬ 
tions is given in lj4l including the single-epoch magnitude 
depths, data quality evaluation, efficiency vs. S/N, and 
the anomalous scatter of flux measurements for objects 
on bright galaxies. The efficiency of multiple detections 
required for a transient is described in 33 including the 
discovery efficiency and the classification efficiency. In 33 
we compare the simulation to data in a preliminary pho¬ 
tometric analysis. Comparisons with SNLS and Diffimg 
limitations are discussed in 33 and we conclude in 33 

2. OVERVIEW OF THE DARK ENERGY SURVEY AND 
TRANSIENT SEARCH 

The Dark Energy Survey includes a wide-area 
5000 deg^ optical survey in the southern celestial 
hemisphere and a dedicated transient search over 
27 de g^, both using the Dark Energy Camera (DE- 
Cam: iFlaugher et al.l[2015ll . DECam is mounted on the 
Blanco 4-m telescope at the Cerro Tololo Inter-American 
Observatory (CTIO) and the data are processed by 
the DES data management system (iSevilla et al.l 120111 : 
IMohr et al.ll2012l : lDesai et al.ll2012f) at the National Cen¬ 
ter for Supercomputing Applications (NCSA). The 570 
Megapixel DECam has a 3 deg^ field of view and is com¬ 
posed of 62 science-image CCDs, each with 2k x 4k pix¬ 
els, and 8 CCDs for guiding. After accounting for CCD 
gaps and two non-functioning CCDs, the active field of 
view is 2.7 deg^. 

The transient search is performed in 10 “SN fields” (27 
deg^) that are repeatedly observed in the g,r,i,z pass- 
bands. We refer to this part of the survey as DES-SN. 
Eight of these fields are observed with few-minute ex¬ 
posure times and are referred to as “shallow” fields; the 
remaining two “deep” helds are observed much longer 
(Table [T]). Defining the AB magnitude-depth as the 
mag where the Diffimg single-epoch detection efficiency 
has fallen to 50%, the shallow and deep field depths are 
^ 23.5 and ~ 24.5, respectively, and the depth in each 
band is the same. The SN portion of the DES observing 
strategy is that the wide-area survey transitions to ob¬ 
serving SN helds when the seeing is above 1.1", and, in 
addition, any SN held (in any band) which has not been 
observed for 7 days is scheduled with the highest observ¬ 
ing priority regardless of the seeing. This 7-day trigger 
typically results in better data quality compared to the 
1.1" trigger. For each SN held, the pointing at a repeat 
visit is the same to within a few arcseconds. Additional 
dithered observations from the wide-area survey are used 
for the inter-calibration of the CCDs. 

On a given night, the number of consecutive exposures 
(A^expose) varies with band and held as shown in Table [TJ 
A^expose = 1 for the shallow g, r, i bands where the sky 
level is well below saturation. In the deep helds (and 
shallow z band), A^expose > 1 to limit the sky level to be 
well below saturation within each exposure. For a shal¬ 
low held, each observing block is scheduled for all four 
bands and takes ^ 20 minutes with overhead. For a deep 
held, observing all exposures in each band takes about 


2 hours, and would be difficult to schedule such a long 
block within the constraints of the global DES observing 
strategy. Each deep-held band is therefore scheduled in¬ 
dependently; the total exposure time per epoch (Texpose) 
is 10 minutes in the g band, and more than an hour in 
the 2 band (Table [T]). 


TABLE 1 

Exposure Summary of DES-SN Fields. 



Band 

Central A 

(A) 

-^expose (sec)^ 

per Epoch 

Template 

(Afepoch)*’ 

Shallow 

9 

4830 

1 X 175 = 175 

8.0 


r 

6430 

1 X 150 = 175 

8.5 


i 

7830 

1 X 200 = 200 

9.3 


z 

9180 

2 X 200 = 400 

9.3 

Deep 

9 

4830 

3 X 200 = 600 

5.5 


r 

6430 

3 X 400 = 1200 

7.5 


i 

7830 

5 X 360 = 1440 

9.3 


z 

9180 

11 X 330 = 3630 

8.3 


^ Yexpose X iexposel — lexpose- 

^ Averaged over fields. The average template exposure time is 

(.^epoch) X ^expose* 

The ten SN helds are divided into four groups of ad¬ 
jacent helds: three C helds that overlap the Chandra 
deep helds, three X helds that overlap the XMM-LSS 
helds, two S helds that overlap SDSS stripe 82, and two 
E helds that overlap the ELAIS SI held. The held loca¬ 
tions were chosen based on (1) visibility from CTIO, (2) 
visibility from telescopes in the northern hemisphere to 
perform follow-up spectroscopy of live targets, (3) galac¬ 
tic extinction, (4) avoiding overlap with extremely bright 
stars, (5) overlap with pre-existing galaxy catalogs and 
calibration. A summary of each held and its location is 
given in Table [3 The maximum nightly data volume 
from observing all ten helds is 170 GB, corresponding to 
just over 5000 CCD images. The average data volume in 
a typical night corresponds to a few helds. In addition to 
the SN held observations, DES-SN makes use of the ex¬ 
tensive calibration data 1 T3.1.1|) taken as part of survey 
operations. 


TABLE 2 

DES-SN FIELD NAMES AND LOCATIONS 


Field 

Deep or 
Shallow 

Field Center (deg): 
R.A. DecL 

Yvisit*" 

g/r/i/0 

Cl 

shallow 

54.2743 

-27.1116 

29/30/30/30 

C2 

shallow 

54.2743 

-29.0884 

28/28/27/28 

C3 

deep 

52.6484 

-28.1000 

25/23/28/27 

XI 

shallow 

34.4757 

-4.9295 

26/27/27/27 

X2 

shallow 

35.6645 

-6.4121 

26/26/25/24 

X3 

deep 

36.4500 

-4.6000 

21/20/22/24 

SI 

shallow 

42.8200 

0.0000 

29/29/28/28 

S2 

shallow 

41.1944 

-0.9884 

27/28/28/28 

El 

shallow 

7.8744 

-43.0096 

27/26/27/26 

E2 

shallow 

9.5000 

-43.9980 

26/26/26/27 


^Number of single-epoch visits to each field in each band, in Yl. 

Science Verihcation (SV) took place 2012 November 
through 2013 January, with the goal of ensuring that 
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the DECam performance meets the DES science require¬ 
ments. During the beginning of SV, the SN fields were 
observed to obtain initial calibrations and to build tem¬ 
plates. The latter part of SV was used to test Difflmg. 
Nominal survey operations began in the Fall of 2013. The 
first season (2013 August to 2014 February) is referred 
to as Yl, and the second season (2014 August to 2015 
February) is referred to as Y2. 


3. THE DIFFERENCE-IMAGING PIPELINE 

All images taken with DECam at CTIO are transferred 
to NCSA and run through the detrending process to pro¬ 
duce images suitable for higher level analyses. For each 
exposure, all CCDs on the focal plane are processed as a 
single unit where bad pixels are masked and corrections 
are applied for bias, flat-held illumination, pupil ghost, 
cr osstalk, linearity, a n d overscan. More d etails are given 
in iMohr et al.l (1201211 : iDesai et Id] (|2012ll and references 
within. The detrending process is virtually identical for 
the SN fields and the wide-area survey, and it is similar 
to a community pipeline used to process DECam data 
for non-DES observers. 

For the SN fields, Difflmg is run after the detrending 
and a schematic overview is shown in Fig. [TJ In contrast 
to detrending, Difflmg is run independently for each 
CCD in order to simplify the distribution of jobs among 
CPUs. Many of our Difflmg sta ges use publicly ay ail- 
able Terapix/AstrOmatic codej^ l)Bertin et al.ll2002tl in- 
cluding SCAMP (lB ertinll2006ll for astrometry, SExtractor 
(iBertin &: ArnoutsI 1 199611 to find objects, PSFEx (iBertinI 
[ml" to determine the position-dependent PSF, and 
SWarp to sum individual exposures (to make “coadds”) 
and to align template images to search images. The sub¬ 
sections below describe Difflmg in more detail. 

3.1. Pre-survey Observations and Analysis 
3.1.1. Calibration 

For the results presented here, the calibratior0 is de¬ 
termined from data taken during SV. The calibration is 
needed to determine magnitudes for detected transients, 
which are used to select transients of appropriate bright¬ 
ness for spectroscopic observations. The calibration is 
also used to convert the fake magnitudes into fluxes in 
CCD counts. 

During nightly operations, DES typically observes a set 
of 3 standard star fields corresponding to low, intermedi¬ 
ate and high airmass. These observations are done dur- 
ing evening twilight , and again during morning twilight 
([Tucker et al.l 1200^ G.Bernstein et al. 2015, in prep). 
These standard star fields are mostly in SDSS stripe 82, 
but supplemented with additional fields, mostly at deck 
~ —45° to —40°. The stars in these fields, which we re¬ 
fer to as secondary standard stars, have had their magni¬ 
tudes transformed into the defined DES “natural” system 
in which the color terms are close to zero. For photomet¬ 
ric nights, these well calibrated secondary standard stars 
are used to determine a nightly calibration consisting of 
zero points (ZPs), atmospheric extinction coefficients. 


and color terms needed to transform the photometry 
from the individual DECam CCDs to the defined DES 
“natural” system. A set of standard stars within each 
of the ten SN fields, referred to as “tertiary standards,” 
were calibrated from data taken during the SV period un¬ 
der photometric conditions, using exposures centered on 
the SN fields plus additional dithered exposures from the 
wide-area survey; this resulted in typically 100-200 well- 
calibrated tertiary standard stars per CCD-area in the 
SN fields (K.Wyatt et al. 2015, in prep). These tertiaries 
were used to calibrate the template images for Difflmg 
f T3.1.2[l . and this calibration is transferred to each tran¬ 
sient magnitude. The relative calibration between DES 
fields over large areas has been checked using the stel¬ 
lar lo cus regression method (iHigh et al.ll2009HKellv et ^ 
l2014fl . where consistency of colors is verified at the 2% 
level. The absolute calibration has been checked at the 
2% level using very short DES exposures on a handful 
of spectrophotometric standards measured by the Hub¬ 
ble Space Telescope. While this early calibration meets 
some of the DES requirements, extensive efforts continue 
to significantly improve the calibration for analysis. 

3.1.2. Templates 

For Yl we constructed deep coadded templates from 
the Y2 season, while the calibration is from SV. Start¬ 
ing with the image that has the lowest sky noise (crgl^Y); 
up to 10 epochs are selected with the smallest PSF that 
have sky noise less than 2.5 ■ The average num¬ 

ber of coadded epochs per band is shown in Table [1] 
along with the total exposure time per epoch. In the 
deep-field i band, for example, the templates include a 
CCD-average of 9.3 epochs which corresponds to a total 
exposure of 3.7 hr. With an average of 8 coadded epochs 
per template, the image-subtracted sky-noise (ctsky) is 
only 6% higher compared to using an ideal template with 
infinite S/nF^ 

Calibrated tertiary standards f T3.1.111 are used to de¬ 
termine the ZP for each exposure, and the pixel flux 
values in each exposure are re-scaled to a common zero 
point, ZP = 31. 1928(3 The coadded templates are com¬ 
bined with a weighted average of each exposure, and the 
weight within each CCD is fixed to the inverse of the 
average sky-variance. 

The astrometric alignment was done in two steps. 
First, the exposures wer e aligned to the USNO-Bl cat¬ 
alog ([Monet et al.ll200§l and then coadded to produce 
an intermediate set of templates in which the alignment 
is good to ^ 100 mas, or 0.4 pixel. These intermediate 
templates were used to produce an internal DES catalog 
based on SExtractor output. Next, the exposures were 
re-aligned to this internal catalog, resulting in ~ 20 mas 
(0.08 pixel) precision. After this final astrometric align¬ 
ment, the exposures are coadded again to produce the 
final set of templates. 


f^SKY/<rSKY (ideal) ~ ^{1 -|- 1/Aftemplate) where Aitemplate is 
the number of coadded templates. 


http://www.astromatic.net 

The global DES calibration plan is available at ZP= 25 + 2.51ogio(300), where 25 is a nominal ZP per second 

http://des-docdb.fnal.gov:8080/cgi-bin/ShowDocument?docid=6584 ^veasidrS©0 sec is a reference exposure time. 
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Fig. 1.— Schematic overview of the SN-field processing: detrending (left panel) and Diffimg (middel+right panels). The thin-lined 
boxes refer to operations on individual exposures; the thick-lined boxes refer to operations on coadds. Astromatic.net codes are shown in 
parentheses. 


3.2. Single-CCD Processing 
3.2.1. Astrometry 

During SV and Yl, the astrometric solution was ob¬ 
tained for the entire focal plane in the detrending pro¬ 
cess, using_the_SCAMP_Qrogr am and the UCAC-4 cata¬ 
log (|Zacha.rias et al.l [201311 as an astrometric reference. 
While this worked well for the wide-area survey, there 
were sometimes very poor solutions in the SN fields lead¬ 
ing to errors up to an arcsecond. We suspected that 
bright saturated objects contributed to this problem be¬ 
cause of the longer exposure times in the SN fields com¬ 
pared to the wide-area survey. 

After Yl, two astrometry updates were incorporated. 
First, we switched to usin g a fainter referen ce catalog 
in the SN fields, USNO-B (|Monet et al.l 1200^ . Second, 
rather than using SCAMP to separately find an astromet¬ 
ric solution for the search and template images, we used 
the SCAMP feature allowing a joint astrometric solution 
for the search and template images. While the absolute 
astrometric precision of the USNO-B catalog (250 mas) 
is worse than UCAC-4 (60 mas), the second change en¬ 
sures good astrometric alignment (< 30 mas) between 
the search and template, which is critical for good sub¬ 
tractions. 


These changes were not incorporated into the detrend¬ 
ing process, and were instead added to Diffimg. Since 
Diffimg is designed for single-CCD processing, a SCAMP 
solution is obtained separately for each CCD rather than 
over the focal plane. The astrometric changes worked 
significantly better, but a few percent of the processed 
CCDs still suffered catastrophic failures in the astromet¬ 
ric solution. As a final refinement to eliminate these 
catastrophic solutions, we used our own DES data to 
construct a reference catalog 1 >13.1.21) . 

3.2.2. Overlaying Fakes onto Images 

In the next stage, two classes of fake point sources 
are overlaid on the CCD image. The first class consists 
of four 20th mag fakes in each band (hereafter called 
“MAG20” fakes) overlaid in random locations away from 
masked regions. The resulting S/N from the Diffimg 
flux measurements is part of the data quality evaluation 
(®. 

The second class of fakes, “SN fakes,” consists of SN la 
light curve fluxes overlaid onto the CCD image near real 
galaxies. The fake SN la light c urve magnitudes ar e 
generated by the SNANA simulation (iKessler et al.ll200^ , 
and include true parent populations of stretch and color , 
a realistic model of intrinsic scatter (|Guv et al.l 120101 : 
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iKessler et alJl201^ . a redshift range from 0.1 to 1.4, and 
a galaxy location chosen randomly with a probability 
proportional to its surface brightness density. All fake 
SN la light curves are generated and stored prior to the 
start of the survey in order to simplify the overlay soft¬ 
ware in Dif f Img. The fake SN la flux added to the image 
is determined by a ZP based on the comparison of cal¬ 
ibration star magnitudes with their fluxes recovered by 
SExtractor. The SN flux is spread over nearby pixels 
using the PSF found by the program PSFEx, and the flux 
in each pixel is smeared by random Poisson noise. 

Ideally, fake SNe would be overlaid onto a duplicate 
set of images so that images with and without fakes can 
be processed separately. For DES-SN we did not prepare 
for this duplication, and therefore care is taken to avoid 
consuming too many galaxies with fake SNe that can 
overlap real transients and cause them to be undetected. 
Fig. [2] shows the fraction of catalog galaxies populated 
by fake SNe la as a function of redshift; the redshift 
distribution has been sculpted to ensure adequate low- 
redshift fakes for monitoring without populating more 
than a few percent of the galaxies at the low and high 
redshift ranges. At a given epoch, the average number 
of overlaid fakes per CCD is ^ 20. Most of the overlaid 
fakes are far from peak or at high redshift, and thus only 
about 1 /3 of these are bright enough to be detected. 

There are a few caveats regarding the selection of 
galaxies and the placement of the fake. First, simulated 
SNe la are matched to a real galaxy based on the galaxy 
photo-z (zphot) since we do not have a sufficiently large 
catalog based on spectroscopic redshifts. To avoid ex¬ 
treme photo-z outliers, we remove galaxies that are ex¬ 
ceedingly bright or faint for its Zphot value by requiring a 
brightness-redshift constraint for both the r and i band 
magnitudes [rrir^i), 

M('^phot) 23 < ^ /^(^phot) 16 : (1) 

where /i(2;phot) is the distance modulus for a flat ACDM 
cosmology with Da = 0.7 and Hq = 70 (km/s)/Mpc. 
This caveat has negligible impact because the fakes are 
overlaid over a wide redshift range and a wide range of 
galaxy mags. 

The second caveat is that the surface brightness pro¬ 
file is assumed to be Gaussian (Sersic index = 0.5) rather 
than a more general sum of Sersic profiles such as a bulge 
plus disk component. This overly simplistic profile re¬ 
sults in fakes placed preferentially near the galaxy cores 
with inadequate sampling of the disk tails. While this 
feature may actually help monitor subtraction problems 
on galaxies, it can result in biased estimates of quantities 
that depend on the distance to the galaxy core, such as 
measuring the fraction of SNe correctly matched to its 
host galaxy. 

The final caveat concerns masking of bad pixels. While 
the placement of fakes is independent of the masking, the 
efficiency analysis presented here ignores fakes in which 
more than 10% of their PSF-weighted pixels are masked; 
7% of the fakes are therefore discarded. For analyses 
requiring the absolute efficiency, such as rates, we can 
impose masking cuts on the data, or perform additional 
fake studies to include the effects of masking. 
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Fig. 2.— As a function of redshift, fraction of catalog galaxies 
in the SN fields with an overlaid fake SN la. 


3.2.3. Image Coadding and Subtraction 

The program SWarp is used to co-add the search ex¬ 
posures, and to remap the template to be aligned with 
the coadded search image. For image subtraction, we 
use a modified version of the difference imaging program 
hotPants. Our version is an attempt to improve the 
performance, and is based on the imn lementatiorF^ of 
A. Be cker, which uses the algorithm of lAlard fc Luntonl 
(|1998l hereafter AL98), and uses some of their original 
code. The basic approach in hotPants is to transform 
one image (which we call a template with pixel values 
tx,y) SO that it can be subtracted pixel by pixel from 
another image taken at a different time and under differ¬ 
ent observing conditions. This linear transformation is 
described by 


y-y =+i~ 

‘C= E 

y-y'^-r 


x—x =+r 

E 

x—x' ——r 


^x{x — x')y(^y—y')^x'y' ( 2 ) 


where t'^-y is the convolved image which is subtracted 
pixel-by-pixel from the unconvolved image. The main 
computation in hotPcUits involves the determination 
of the values of the kernel of the transformation 
kx{x-x')y{y-y’)- The parameter r is the size of the kernel 
and X, y, x' and y' are the pixel coordinates. In general, 
one should add a constant term to Eq[5J but our version 
makes a global background subtraction of the images be¬ 
fore determining the kernel. 

The kernel is assumed to vary slowly over the image 
and this variation is described by a polynomial; 

h , 

’^x(x-x')y(y-y') — i^(x-x'){y-y') 

+ XK^x-x'){y-y') 

+ y^ll-x'){y-y’) + ( 3 ) 

The AL98 algorithm allows a polynomial of arbitrary 
order, but since we process each CCD separately our 
hotPants version includes only linear terms. 

A major difference between our version and the 
AL98 algorithm lies in the parameterization of the 
kV ,w AL98 parameterize the kernel as an arbi- 
trary number of Gaussian functions of fixed width mul¬ 
tiplied by polynomials whose coefficients are parameters 


http://www.astro.Washington.edu/users/becker/v2.0/ 
hotpants.html 









to be fit. In routine use, the number of polynomial coeffi¬ 
cients is large and comparable to the number of pixels in 
the ker nel. Inste a d, we have chosen a method similar to 
that in iBramici] (j2008[l in which the pixel values in the 
core of the kernel are fitted without the use of a function 
to parameterize them. We have, however, retained the 
Gaussian function for the pixels at the edges of the ker¬ 
nel: the Gaussian form is useful for cases where a large 
kernel is needed to match images with very poor seeing. 
While our approach seems more transparent in terms of 
understanding the fit parameters, we do not have solid 
evidence that our parameterization results in better sub¬ 
tracted images. 

3.2.4. Detections and Candidates 

We first measure the PSF to define the detection profile 
we are searching for, and then PSF-like objects on the 
subtracted image are found by SExtractor. Selection 
requirements in Table [3] are applied to reduce the number 
of artifacts. An object satisfying these requirements is 
referred to as a “detection.” 

A “raw candidate” is defined when two or more de¬ 
tections have measured positions matching to within 1". 
The two detections can be in the same band or differ¬ 
ent bands, or on the same night or different nights. All 
raw candidates are saved, which includes moving objects 
such as asteroids and KBOs. Requiring detections on 
separate nights 1 H3.3I) is used to reject moving objects. 

TABLE 3 

Selection Requirements on SExtractor Detections in a 
Subtracted Image 


(1) S/N > 3.5, although the effective S/N cut from 
SExtractor is higher (~ 5) as shown in Fig. [8] 

in 35 X 35 pixel stamp around the detected object: 

(2) fewer than 200 pixels with a flux less than —2a below zero, 

(3) fewer than 20 pixels with flux less than —4a below zero, 

(4) fewer than 2 pixels with flux less than —6a below zero. 

(5) detection not near object in veto catalog containing 
80,000 stars with r-band mag < 21. 

Veto radius is mag-dependent, and total vetoed area 
over all 10 fields is 0.63 deg^, or 2.4% of the area. 

(6) for co-added images, cosmic ray rejection based on 
consistency of detected object on each exposure. 

(7) detected object profile is PSF-like based on the 

SExtractor SPREAD_M0DEL variable f ^esa^^a^^f^^l _ 

(8) SExtractor A_IMAGE < 1.5 X PSF 


3.3. Post-processing 

In addition to the single-GGD operations, there are 
post-processing steps that operate on all fields and 
GGDs, and continually update the candidate properties. 
A few percent of the events land on a GGD in two over¬ 
lapping helds, and thus single-GGD processing is not a 
useful concept when constructing candidates from mul¬ 
tiple observations. 

In some past surveys, as well as the start of DES-SN, 
the first post-processing step was to perform a visual in¬ 
spection of each detection in order to reject subtraction 
artifacts that produce false detections. In DES-SN we 


use a new machine learning based code to replace hu¬ 
man scanning; this “autoScan” program is described in 
detail in GI5. The algorithm makes use of the supervised 
machine learning technique Random Forest. The train¬ 
ing sample includes nearly 900,000 Diffimg detections, 
half of which were flagged as artifacts by human scanners 
and the other half are detections of fakes. For each de¬ 
tection, the inputs to autoScan include a 51 x 51 pixel^ 
detection-centered stamp from the search, template, and 
subtracted images. The flux and uncertainty from each 
pixel on these three stamps contributes ^ 15, 000 pieces 
of information. However, rather than using the pixel- 
level information we found that autoScan performs bet¬ 
ter and faster using 37 high-level features computed from 
the stamps. The three most important features are (1) 
ratio of PSF-fitted flux to aperture flux on the template 
image, (2) mag-difference between the detection and the 
nearest catalog source, and (3) the SPREADJfODEL output 
from SExtractor. 

For each object, the autoScan program returns a score 
between 0 and 1, where 0 corresponds to an obvious 
artifact and 1 is for a high-quality detection. While 
autoScan could have been applied before making raw 
candidates, we have so far been conservative and apply 
the autoScan requirement here in the post-processing in 
order to fully monitor the autoScan performance. 

The first post-processing step is to define “science can¬ 
didates,” a detection on two distinct nights, each satis¬ 
fying the autoScan requirement. Science candidates are 
the official product of Diffimg, and as more epochs are 
acquired these candidates are repeatedly analyzed to se¬ 
lect targets (object and host galaxy) for spectroscopic 
observations. If there is a future science case requiring 
single-night detections, we can recover the raw single¬ 
night candidates; the caveat is that during survey oper¬ 
ations, only the 2-night science candidates are selected 
for spectroscopic observations. 

The next post-processing stage is to match each science 
candidate to a host galaxy, which is later targeted for 
a spectroscopic redshift. We use th e “direction a l-light - 
radius” (cIlr) method described in ISako et ffil (I20I4I1 . 
Gurrently the galaxy profiles are approximated by a 
Gaussian (Sersic index = 0.5), and will eventually be 
updated with profile fits to an arbitrary Sersic index. If 
there are multiple nearby galaxies within 4 x duR they 
are all flagged to acquire a spectroscopic redshift. 

The next post-processing stage, “forced photometry,” 
computes the PSF-fitted flux and its uncertainty for each 
observation since the start of the observing season, re¬ 
gardless of whether there was a detection. The flux and 
uncertainty are computed at the same coordinates (R.A., 
deck) on each subtracted image, and the coordinates 
are computed as the weighted average from each detec¬ 
tion. This stage allows recovering small fluxes just be¬ 
low detection threshold, and fluxes consistent with zero, 
in order to construct complete light curves. Ideally the 
autoScan program would be used to flag bad subtrac¬ 
tions that could lead to badly measured fluxes. How¬ 
ever, while the autoScan results exist for detections, we 
do not have the infrastructure to run autoScan on non¬ 
detections in a manner analogous to the forced photom¬ 
etry. In addition, autoScan would need additional train¬ 
ing to accept subtractions with no significant detection. 
As an alternative to autoScan, forced photometry mea- 
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surements are rejected from light curve fitting (below) if 
(1) the PSF-fitted flux and aperture flux differ by more 
than 5(7, or (2) within a radius there are 2 or more 
pixels with S/N < —6. 

The final po s t-processing stag e is to use the SNANA pro¬ 
gram PSNIIO (|Sako et al.ll20li() to perform photometric 
classification by comparing each candidate light curve to 
a series of photometric griz light curve templates con¬ 
structed on a redshift grid for (1) SN la, (2) CC type 
II, and (3) CC type Ib/Ic. For each candidate-template 
calculation, we discard up to two epochs with the 
largest contribution (if above 10). This outlier rejec¬ 
tion helps to avoid bad fits from a few poorly measured 
forced-photometry fluxes, particularly on bright galaxies 
as described in i]ll A relative probability is computed 
from each x^; a-nd a Bayesian probability is computed 
for each SN type; the largest probability (Pmax) deter¬ 
mines the type and redshift. If Pmax < 0.5, or the best 
fit x^ is poor, the candidate type is flagged as unknown. 
The probability for each type and the estimate of peak 
magnitude contribute to the spectroscopic target selec¬ 
tion process (' il3.4.2|) . 

3.4. Spectroscopic Target Selection 

While spectroscopic target selection is outside the 
scope of Dif f Img, here we give a brief description to give 
a more complete picture of the DES-SN program. The 
two components of spectroscopic targets, host redshifts 
and live transients, are described below. 

3.4.1. Host Galaxy Redshifts 

The large numbers and faint magnitudes of SNe dis¬ 
covered in DES-SN overwhelm the available resources 
for spectroscopically classifying each candidate. How¬ 
ever, we can efficiently use multi-fiber spectroscopic 
resources to measure an accurate host-galaxy redshift 
for the majority of our SN candidates. Using the 
Anglo-Australian Telescope (AAT), the OzDES program 
(jYuan et al.l I20I5H is a lOO-night spectroscopic survey 
with the 400 fiber Two Degree Field (2dF) instrument 
feeding the dual-beam AAOmega spectrograph. The 
overlap between the field of view of DECam and 2dF is 
nearly complete. With repeat visits to the same source, 
spectra are coadded to enable redshift measurements for 
much fainter galaxies than would naively be expected 
from a 4 m class telescope; redshifts are obtained for 
about half of the 24th mag galaxies {r band). In addi¬ 
tion to targeting host galaxies for SN candidates, OzDES 
also targets a variety of DES sources such as AON to 
derive reverberation mapped black-hole masses, galaxies 
for DES photo-z calibration, white dwarfs for calibration, 
and live transients for spectroscopic typing. 

3.4.2. Spectroscopic Identification of Live Targets 

The spectroscopic selection for live transients is pri¬ 
marily focused on SN la. The selection is based on 
a visual examination of light curves along with PSNID 
probabilities. The phase estimate is used to give higher 
priority to candidates near peak brightness. Highest pri¬ 
ority is given to candidates with peak r band magnitude 
J’peak <20.5 mag (mag-limited) and to candidates with a 

PSNID—Photometric SN Identification 


photometric redshift below 0.2 (volume limited). These 
two samples have large overlap, and are expected to be 
very nearly complete. Lower priority is given to can¬ 
didates over the full redshift range where we expect to 
acquire a spectroscopic typing for ~ 10% of the SN la 
sample. Starting in Y2, transient activity in multiple 
seasons is used to reject AGN-like candidates. 

Telescopes used to spectroscopically confirm tran¬ 
sients discovered by D iff Img include the 3.9-m Anglo- 
Australian Telescope (AAT) at Siding Springs Observa¬ 
tory in Australia, the 8.2-m Very Large Telescope (VLT) 
on Cerro Paranal in Chile, the 9.2-m South African Large 
Telescope (SALT) near Sutherland in South Africa, the 
10.4-m Gran Telescopic Canarias (GTC) in La Palma, 
the Keck 10-m on Mauna Kea in Hawaii, the 6.5-m Mag¬ 
ellan Telescope at Las Campanas Observatory in Chile, 
the 6.5-m MMT on Mount Hopkins in Arizona, the 3- 
m Shane telescope at Lick Observatory in California, 
the 4.1-m Southern Astrophysical Research (SOAR) tele¬ 
scope at Cerro Pachon in Chile, the 8.1-m Gemini-South 
telescope at Cerro Pachon in Chile, and the 9.2-m Hobby- 
Eberly Telescope at the McDonald Observatory in Texas. 

3.5. Statistics Summary 

A summary of the first-season (Yl) statistics for single¬ 
epoch detections is shown in Table |4l The average num¬ 
ber of objects per field found by SExtractor increases 
with the passband central wavelength. In the shallow 
fields there are ~ 100, 000 per field in the g band, in¬ 
creasing to ^ 170,000 in the z band. In the deep fields 
there are 130,000 in the g band, increasing to 270,000 
in the z band. A visual scanning assessment shows that 
more than 90% of these detections are subtraction arti¬ 
facts. Following the SExtractor detections on the sub¬ 
tracted image, there is a significant reduction from the 
selection cuts and autoScan. The selection cuts reduce 
the number of detections by a factor of 3-4 in the g band, 
and a factor of ^ 2 in the z band. The automated scan¬ 
ning provides a further reduction of a factor of ^ 4 in 
the g band, increasing to an order of magnitude in the 
z band. After all selection requirements and automated 
scanning, the average number of objects per field in Yl 
is ^ 10^ in both the deep and shallow fields, and the 
artifact fraction is ^ 25% as determined from a visual 
scanning assessment. 

To determine the average number of detections per 
square degree for a single-epoch visit (h-detect), the num¬ 
ber of Yl detections (autoScan row in Table |3]) is divided 
by 2.7 deg^ and Wisit from Tabled hdetect « 110 in the 
g band and ps 150 in the z band. 

The total number of raw candidates in Yl, which re¬ 
quires two SExtractor detections passing the selection 
cuts in Table m is 1.2 x 10^. Requiring two detections 
on different nights reduces this slightly to 1.0 x 10®. Re¬ 
quiring the two separate-night detections to satisfy the 
automated scanning reduces the number of candidates to 
7489, or a factor of 13 reduction. Table[5]shows the aver¬ 
age number of candidates per deep field and per shallow 
field. 

Following SExtractor detections, the selection cuts 
and autoScan have a dramatic effect on reducing the 
number of detections and candidates. This is because 
the vast majority of the SExtractor detections are false 
positives, or artifacts of the image subtraction. These 
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TABLE 4 

Number of Non-fake Single-Epoch Detections in the Y1 
Season per 3 deg^ Field (thousands) 


Fields 

Detection 

Stage 

Number of Detections 
(xlO^) 

g r i z 

Deep 

SExtractor 

133 

166 

277 

270 


+ selection cuts^ 

32 

81 

172 

167 


+ autoSccin^ 

8 

8 

9 

12 


autoScan/cuts ratio 

0.25 

0.10 

0.06 

0.07 

Shallow 

SExtractor 

98 

103 

126 

173 


+ selection cuts^ 

29 

26 

55 

92 


+ autoSccin'^ 

8 

7 

9 

10 


autoSccin/cuts ratio 

0.28 

0.27 

0.18 

0.11 


^Includes selection cuts in Table |3] 

'^Includes cuts and automated scanning requirement (G15). 


TABLE 5 

Average Number of Non-fake Y1-candidates per 
3 deg^ Field. 


Candidate 

Acand per Field 

Selection 

DEEP 

SHALLOW 

2 detections (raw cand) 

18830 

10410 

2 nights (without autoScan) 

17460 

8230 

2 nights + autoScan (science cand) 

1040 

680 


artifacts come from a variety of sources, including bright 
stars and galaxies, defective pixels, edges of masked re¬ 
gions, CCD edges, and cosmic rays. Some of these arti¬ 
facts are illustrated in Fig. 1 of G15. The large rejection 
by autoScan costs only a 1.0% loss of fake SNe la candi¬ 
dates, mainly for fakes with low S/N at peak brightness. 
We are therefore confident that autoScan is highly effi¬ 
cient for real astrophysical transients. 

Subtraction artifacts are illustrated in Fig.|3]for a deep 
field image processed by Difflmg. SExtractor detec¬ 
tions failing selection cuts (dashed red boxes) are the 
most clearly evident upon visual inspection, while those 
failing autoSccin (solid red boxes) are more subtle. In 
this example, most of the artifacts are around a few 
bright objects even though most of the bright sources are 
cleanly subtracted. On average, artifacts are ~ 1 mag 
fainter than real transients. To get an estimate of the 
artifact rate for bright sources, ^ 3% of bright fakes 
(mag< 20) fail the detection and autoScan requirements. 
The origin of these artifacts is not understood. 

3.6. Classification Summary 

Here we show the breakdown of PSNID classihcations 
for science candidates f il3.3l) . To avoid the noisiest light 
curves we consider the subset in which three bands each 
have an observation with S/N> 5; this subset is roughly 
half of all candidates. Applying PSNID to the entire light 
curves for the full Y1 sample results in nearly equal clas¬ 
sification fractions (~ 1/3) for SN la, SN CC (mostly 
Type II) and unknown. 

While a full Y1 analysis is relevant after the survey, 
during survey operations PSNID is run on newly discov¬ 
ered light curves that have only a few epochs. To il¬ 


lustrate the real-time PSNID performance, Fig. 0] shows 
the classification fractions as a function of time the light 
curve has been observed. MJDcand is the time when the 
second epoch is detected, or when the object became a 
science candidate. MJDpef represents the current MJD, 
which we take to be 56,600 in this example. The fits in¬ 
clude observations between MJDcand — 20 and MJDref. 
When only the early part of the light curve is available 
for fitting (—5 days in Fig. |4]), about 70% of the candi¬ 
dates are classified as SN la, fewer than 10% as SN CC, 
and the rest are unknown. When fitting 2 months of 
the light curve, more than half of the classifications are 
SN CC. 


3.7. Data Reprocessing 

During Yl, the monitoring of fakes showed a significant 
inefficiency that was traced to severe astrometry prob¬ 
lems as described in il3.2.1l This problem was fixed after 
Yl, and before the start of Y2 operations all of Yl was 
reprocessed in order to recover hundreds of host-galaxy 
spectroscopic targets that had been missed during Yl. 
During Y2, the monitoring of fakes showed good Difflmg 
performance in the shallow fields, but there were still sig- 
nihcant flux-outliers in the deep fields. This problem was 
eventually traced to the program which determines the 
PSF used for calculating PSF-fitted fluxes, and it was 
fixed after Y2. 

Both Yl and Y2 have been fully reprocessed in all ten 
SN fields, with all Difflmg fixes. Results presented in 
this paper are based on the Yl season, using templates 
constructed from Y2 images. The reprocessed results are 
used to discover transients missed during the survey, to 
update the photometric classification with PSNID, and to 
update the host-galaxy target list for measuring spec¬ 
troscopic redshifts. While transients discovered in the 
reprocessing have become too faint to target for spectro¬ 
scopic observations, this is not a serious issue because we 
target only a small fraction of the transients anyway. 

Another subtle change in the reprocessing campaign 
was to fully analyze each exposure in the deep field se¬ 
quences (in addition to the coadd) to improve the KBO 
search. In particular, this reprocessing led to the discov¬ 
ery of one of the two Neptune Trojans in iGerdes et al.l 
(EoH, as well as improved orbital fits for both ob¬ 
jects. We are currently upgrading Difflmg to overlay 
fake KBOs onto the images; these fake KBOs will allow 
measuring the search efficiency, and they will be used to 
develop improved KBO-finding algorithms. 

We do not expect more Difflmg improvements during 
the remainder of DES, unless our monitoring uncovers 
new problems or we improve the subtraction problem 
on bright galaxies as described in il4.31 Even without 
software changes, we may reprocess the data in the future 
using better templates and lower detection thresholds in 
order to improve the depth of the search. 

3.8. Difflmg Processing Time 

Using the IBM iDataPlex Carver computational sys¬ 
tem at NERSCS, we give the processing time for the 

National Energy Research Scientific Computing Center 


















11 



Fig. 3. From DiffImg, a co-added search image (top) and subtracted image (bottom) from a typical night (2013 October 13) in deep 
field C3 for i band. The image size is roughly 2.7' x 4.6^ or about 1/13 the area viewed by a single CCD. Non-fake SExtractor detections 
on the subtracted image (bottom) are highlighted in both images: dashed red boxes for objects failing the selection cuts in Table [3] solid 
red boxes for objects passing these cuts and failing autoScan, and yellow circles for objects passing cuts and autoScan, which are used to 
make science candidates. To set the scale, the brightest masked star has mag m = 11.4; the other masked star has mag m = 15.0. To see 
detections in more detail, Fig. 1 in G15 shows a collection of 51 x 51 pixel^ stamps for search-|-template-|-subtracted images, each centered 
on a detection. 


Diff Img steps in the middle panel of Fig. [T] For a shal¬ 
low field with a single exposure, the processing time for 
a single CCD is ^ 10 minutes, half of which is spent on 
the hotPants program. In the deep fields we perform 
the hotPants subtraction for each exposure as well as 
the coadded image, and thus the processing time scales 
roughly with the number of exposures. For a deep-held 
sequence with 11 z band exposures, the processing time 
for a single CCD is ^ 90 minutes. 

The post-processing steps (right panel in Fig. [T]) run 


serially, and the processing time depends on how long 
the survey has been running. Near the start of a survey 
season the post-processing takes a few minutes, but near 
the end of the season it takes several hours. 


4. Dlfflmg MONITORING-I: SINGLE EPOGHS 

Here we describe monitoring of the single-epoch detec¬ 
tion efficiency and data quality, using both the MAG20 
fakes and the SN fakes processed by Diff Img. 
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Fig. 4.— PSNID classification fraction for SN la and SN CC 
vs. time that the light curve has been observed, for candi dates in 
which three bands have an observation with S/N> 5. See N3.6I for 
explanation of MJDcand MJDref- Zero on the horizontal axis 
corresponds to newer candidates used in the PSNID fits; —50 days 
corresponds to older candidates whose second detection occurred 
50 days earlier and thus have longer light curve coverage in the 
PSNID fits. “Unknown” corresponds to light curves for which PSNID 
cannot determine a SN type. 

4.1. Data Quality Assessment 
The measured S/N from the MAG20 fakes is part of the 
data quality evaluation (See Fig. [5]). We define S/Njjj^g 2 o 
to be the average S/N among all of the (4 x 60 = 240) 
MAG20 fakes overlaid on each exposure, where each S/N 
is the ratio of the PSF-fitted flux to its uncertainty. If 
S/Nmag 2 o < 20 in the shallow fields, or < 80 in the deep 
fields, the exposures are flagged to be retaken. In ad¬ 
dition, an exposure is retaken if the i band PSF width 
(FWHM) at zenith is > 2"; this seeing value is computed 
by correcting the measured PSF for airmass and wave¬ 
length. These criteria for retaking an exposure are a com¬ 
promise between data quality in the SN fields and lost 
observing in the wide-area survey. The largest S/Njjjg^g 2 o 
values are from high-quality data triggered because there 
were no observations within the past 7 days. The lower 
S/Nj^^g 2 o values are typically from data triggered by see¬ 
ing > I.l" and from observations at larger airmass. 

While S/Nj^^g 2 o determine 

if an exposure sequence needs to be retaken, the SN la 
fakes are used to determine complementary information 
about the data quality. For a given epoch, the fakes are 
used to determine the magnitude depth, TOeff=i/ 2 i de¬ 
fined as the mag where the Diff Img detection efficiency 
has fallen to 50%. Figure [6] illustrates the determina¬ 
tion of TOeff=i/2- The ?7iefi=i/2 distribution is shown in 
Fig.[7]for each band, and for deep and shallow fields. The 
variation in meff=i /2 is from the variation in observing 
conditions. 

4.2. Detection Efficiency vs. S/N 
The detection efficiency as a function of S/N (cs/n) 
is a crucial input to the MC simulation ([J5|) and also 
provides another monitoring metric. We do not attempt 
a first-principles calculation of es/N; primarily because 
of the complicated behavior of SExtractor that largely 
defines the detection threshold. Therefore es/N is em¬ 
pirically measured from the fakes as illustrated in Fig |8] 


for the i band. The effective S/N threshold, defined for 
es/N = 0.5, is about 5 in each band and is the same in 
both the deep and shallow fields. Each sub-panel shows 
the nominal es/N curve computed from all of the fake 
data, along with a systematic test based on splitting the 
data into two equal-size samples. The probability of de¬ 
tecting a transient depends on the ZP, PSF, and sky- 
noise through their effect on the S/N, and we expect the 
detection efficiency to depend primarily on S/N. Fig [8] 
shows that there is no unexpected dependence, which 
is important because not all of the selection criteria are 
based on S/N. 

4.3. Anomalous Subtractions on Bright Galaxies 

The final issue is the reliability of forced-photometry 
flux measurements that are used to classify light curves, 
both visually and with fitting programs. The average 
fake fluxes are recovered to within few percent of their 
true values, which is adequate precision since it is smaller 
than the model errors used in light curve fitting. We 
have also checked the reliability of the flux uncertainties, 
and found that these uncertainties are underestimated 
in proportion to the local galaxy surface brightness (SB) 
under the SN location; we refer to this effect as the “SB 
anomaly.” The excess flux scatter can cause problems 
with monitoring and light curve fitting, and thus we have 
modeled this effect in both simulations and fitting pro¬ 
grams m- 

To define the SB, we first sum the template flux at the 
candidate location, using an aperture with 1.3" radius, 
which contains most of the flux for a typical PSF. The 
SB flux is defined as the average flux per square arcsec- 
ond, and the SB-mag (mse) is the corresponding mag¬ 
nitude per square arcsecond. For fakes we characterize 
the quality of the uncertainties using the rms of AF/ap 
(RMSa), where AF is the difference between the mea¬ 
sured (forced photometry) flux and the true flux of the 
fake, and ap is the uncertainty on the forced-photometry 
measurement. Ideally RMSa = 1 in all cases, but we find 
that RMSa increases with SB as shown in Fig. IH] for the 
deep fields and in Fig. [TU] for the shallow fields. For low 
SB (msB > 24), RMSa is very close to unity as expected. 
For the brightest galaxies where msB ~ 20, RMSa ~ 5 
in the deep fields and ~ 3 in the shallow fields. 

Figures [5] and [TO] also show rms vs. msB separately 
for dim fakes with m > 26 (red curve) and for brighter 
fakes with m < 24 (blue curve). The consistency shows 
that this effect depends mainly on the brightness of the 
galaxy and not the transient source. 

5. Diffimg MONITORING-II: SCIENCE CANDIDATES 

While monitoring the single-epoch detection efficiency 
and data quality are important on a nightly basis O, 
the science prospects ultimately depend on the Diffimg 
candidate efficiency and our ability to select spectro¬ 
scopic targets based on a small number of epochs. Here 
we describe the Diffimg monitoring of science candidates 
using SN fakes combined with MG simulations. The ba¬ 
sic idea is to use the MG simulation to predict the SNe la 
efficiency versus redshift, and compare with the true ef¬ 
ficiency measured from the fakes. 

There are two different efficiencies to monitor as a 
function of redshift. The first efficiency is the fraction 
of fakes that become a science candidate (Scand)- As 
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Fig. 5.— For each set of exposures in Yl, the CCD-averaged 
PSF width (i-band at zenith, FWHM, arcsec) is plotted against the 
CCD-averaged S/N from the MAG20 fakes. Exposure sequences 
with points outside the dashed box are retaken, typically the fol¬ 
lowing night. Left panels are for the deep-fields and the right panels 
are for the shallow fields. 




20 22 24 26 28 20 22 24 26 28 

true mag true mag 

Fig. 6.— True mag distribution for fakes from a single epoch. 
Left panel is for deep field C3-r; right panel is for shallow field 
Cl-r. Shaded overlay is for fakes satisfying Diff Img detection re¬ 
quirements. The dashed vertical line shows 3-® defined in 

the text. 

long as fcand IS Optimal, then even if Diff Img measures 
fluxes with many catastrophic outliers an improved of¬ 
fline photometry analysis can make all of the discovered 
light curves useful for science analysis. However, if there 
are too many flux outliers then real-time photometric 
classification becomes more difficult, which complicates 
the selection of spectroscopic targets. 

It is therefore important to monitor a second Diff Img 
efficiency, the fraction of fakes passing the photomet¬ 
ric analysis (fpsNro) used for spectr oscopic targeting , 
which is based on the PSNID program (iSako et al.ll201lll . 
The key component of the PSNID analysis (Table [51) is 
a requirement on the fit probability computed from the 
template-fit therefore even a few measured fluxes 

that are highly discrepant from their true values can 
cause PSNID to reject the light curve. Up to two highly 
discrepant fluxes (w.r.t. the fit) are rejected, allowing 
for a small level of subtraction problems. In summary, 
simply discovering an event is not adequate unless the 
flux measurements are of sufficient quality to perform 


Fig. 7.— Distribution of Tneff_x /2 each passband, determined 
with fakes. Each entry is from one epoch. Left panels are for the 
deep fields; right for the shallow fields. 



Fig. 8. — Single-epoch detection efficiency (cs/n) vs. S/N, as 
measured with fakes. The solid-filled circles are computed from 
all the data, and is the same in each panel. The solid and dashed 
curves correspond to splitting the sample into roughly two equal 
bins for zero point (top), PSF (middle) and sky noise (bottom). 

light-curve template fitting without suffering significant 
inefficiency. 

Details of the MC simulation are given in Appendix K1 
and here we give a brief overview. The MC simnlation 
uses the observed cadence, and the simulated flux and 
noise are computed from the observing conditions at each 
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Fig. 9.— For the 2 deep fields in each pass band, RMS of AF/ap 
as a function of the galaxy surface-brightness mag (mgs) for fakes 
as defined in the text. AF is the difference between the true and 
measured fake flux, and ap is the uncertainty. The horizontal 
dashed line through 1 shows the expected value if Diff Img correctly 
determines the flux uncertainties. The dotted-red curve is for very 
faint fakes with SN mag m > 26; the dashed-blue curve is for fakes 
with SN mag m < 24. The model for this effect depends on 
mgB and not on the transient mag. 
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Fig. 10.— Same as Fig|9] but for the 8 shallow fields. 

epoch: ZP, PSF, sky noise, CCD gain. 

While the cadence information is trivially obtained 
from survey observations, the MC simulation also needs 
two inputs based on the fakes processed by Difflmg. 
First, we use the efficiency vs. S/N (c^n) measured in 
each passband, and illustrated in Fig. |3 for the i band. 
Since there is good agreement between the deep and shal¬ 
low fields, we use the same cs/n function in all fields. 

The second input from the fakes is a model for the SB 
anomaly, the anomalous flux uncertainty that increases 
with the local surface brightness. The galaxy Sersic pro¬ 
file in the simulation is used to analytically compute 
msB, and the RMSA-versus-msB curves in Figures [9] and 
[TOl are used to scale the sky noise as a function of pass- 
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TABLE 6 

PSNID Analysis Requirements for fpsNiD Study 


Category 

Requirement 

Sampling 

5 or more observations. 

3 bands with at least one S/N> 5 observation. 

An observation with Tobs < days.^ 

An observation with Tobs > +5 days. 

Fit-x^ 

Fit prob Pfit >0.1^ 

Reject up to two 3.16cr data-fit outliers (Ay^ > 10). 

Typing 

Best fit template (among la, II, Ib, Ic) is Type la. 


'‘Tobs is the observer-frame time since the epoch of peak bright¬ 
ness. 

Pfit is calculated from x^/dof. Because of the large PSNID model 
errors, the true chance of finding Pfn < 10% for SNe la is ~ 1%. 

band, and as a function of deep or shallow field. These 
same RMSA-versus-msB curves are used in the PSNID 
analysis to scale the flux uncertainties. The PSNID anal¬ 
ysis results in 7662 fakes passing the selection criteria in 
Table [6] (includes all 10 fields), and a similar number of 
SNe la from the MC simulation. 

FigurefTTlshows the science-candidate efficiency (fcand) 
and PSNID-analysis efficiency (fpsNio) as a function of 
redshift for one shallow field in each group. The analo¬ 
gous deep-held plots are shown in Fig. [121 In the shallow 
helds, fcand — 1 for redshifts z < 0.5, and falls to 50% 
at z ~ 0.7. In the deep helds, fcand — 1 for redshifts 
z < 0.8, and falls to 50% at z cs 1.1. The overall agree¬ 
ment is good between the fakes and the MC simulation. 
While we might have expected the SB anomaly to affect 
the discovery of lower redshift SNe that preferentialy lie 
on brighter galaxies, we hnd that the low-redshift effi¬ 
ciencies are ~ 100% and thus the SB anomaly has a neg¬ 
ligible impact on discovering SNe la. The SB anomaly 
and its impact are discussed further in 117.21 The most 
notable discrepancy is in fcand for redshifts z > 1.2 in the 
C3 deep held, and SpsNiD for redshifts z > 0.8 in both 
of the deep helds (Fig. IT^ . Finally, it is worth noting 
that prior to the hnal reprocessing the fake efficiencies 
were signihcantly worse than the MC prediction for the 
reasons described in 113.71 

5.1. What are the Science Candidates ? 

Here we give a very approximate breakdown for the 
7500 science candidates discovered by Difflmg in Yl, 
where each candidate requires a Difflmg detection on 
2 separate nights with no other selection requirements. 
First we use our MC simulation to predict the SN contri¬ 
bution (la-bCC; see Appendix and we include events 
that reach peak brightness well before and after the Yl 
season. We hnd 2000 ±300 SNe, where the uncertainty is 
from the rate measurements, and nearly 60% of the SNe 
are Type la. This SN contribution corresponds to about 
27% of the candidates. 

A non-astrophysical candidate, or artifact, is dehned 
as a candidate in which more than half of the detections 
fail the automated scanning requirement (IJSTSjand G15). 
Using this arbitrary but illustrative dehnition, ~ 30% 
of the science candidates are artifacts (i.e., ~ 2300 in 
Yl), compared with 1.5% of the fakes. These artifacts 
become a science candidate because of the relatively loose 
requirement of only 2 detections passing the selection 
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requirements and automated scanning. The relatively 
small number of artifacts does not cause problems during 
survey operations, and thus we choose to reject them 
with offline analysis software rather than trying to reduce 
the number of science candidates. 

For the remaining science candidates, a preliminary 
assessment of the OzDES spectral classifications shows 
that they are mostly AGN and variable stars. 
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Fig. 11. — For one shallow field in each group, the Y1 efficiency 
vs. redshift is shown for fakes processed by Diffimg (black dots), 
and the MC prediction (red histogram). Left panel is the efficiency 
for becoming a science candidate (fcand); right panel is the PSNID- 
analysis efficiency (fpsNlD) defined in Table |6] 


6. REALITY CHECK: DATA-MC COMPARISON 

Since the Diffimg results presented so far are based 
on fakes and simulations, here we perform a reality check 
and compare the SNANA-based MC simulation to Y1 data, 
where the MC simulation is a mix of SNe la and CC SNe 
as described in Appendix Recall that the MC simu¬ 
lation has input from fakes processed by Diffimg, but 
there is no tuning with real science candidates. Here we 
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Fig. 12.— Same as Fig. 1111 but for the two deep fields. 


make a data-MC comparison for the photometric redshift 
distribution (zphot) of a photometrically selected SN la 
sample, using only the SN light curve information. We 
do not use any spectroscopically confirmed typing infor¬ 
mation, nor do we use any host-galaxy redshifts. 

For this comparison we do not use the PSNID selec¬ 
tion criteria in Table [H Instead, we use a more strin¬ 
gent analysis designed to photometrically select a highly 
pure SN la sample. We fit both the data and MC sam¬ 
ples with the SALT-II model (iGuv et al.ll201(]D usiii g 
the photo-z technique described in iKessler et al.l (|2010bf) . 
Finally, the SALT-II fit parameters are used in a near- 
est neighbor (NN ) analysis similar to that described in 
iSako et ^ ( 20I4D . Details of the analysis are given in 
Appendix [B] and the resulting Zphot comparison is shown 
in Fig. [131 The high-redshift roll-off in the Zphot distri¬ 
butions is mainly from the requirement that three bands 
each have an observation with S/N> 5. The overall 
agreement is reasonable, except for z > 1 in the deep 
fields. This data-MC discrepancy will be monitored as 
we continue to improve photometric classification meth¬ 
ods and the simulation. 
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Fig. 13.— ^phot distribution in Y1 from a photometric analysis 
(Appendix [BJ. Left panel is for the 2 deep fields; right is for the 
8 shallow fields. Black points are the data; red histogram is the 
MC prediction, re-scaled to have the same number of entries as the 
data. 
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7. DISCUSSION 

7.1. Comparison of Search with SNLS 

Here we make some rough performance comparisons 
between the SNLS and DES-SN deep field search for 
SNe la. These two surveys have similar depths and pass- 
bands, and they each measured their efficiency with fake 
SNe la overlaid on images. While the DES-SN trigger re¬ 
quires 2 epochs in any band, the SNLS trigger requires a 
single detection in the Megacam band. For the single¬ 
epoch detection efficiency. Fig. 9 of PIO shows that the 
magnitude at 50% efficiency is mefr=i /2 = 24.3 in the 
band for an exposure time of 3640 sB This depth is very 
similar to our average DES i band depth, mefr=i /2 = 24.5 
(Fig. [7]), using 1440 s exposures. 

PIO also measure the efficiency vs. redshift for finding 
fake SNe la. Both the PIO and SNANA simulations predict 
the observed color and stretch distribution for SNLS, and 
thus the two simulations are consistent in describing the 
parent populations of stretch and color. Fig. 10 of PIO 
shows that fcand = 50% a.t z ce 0.95, slightly below the 
corresponding DES-SN redshift z ~ 1.1. 

7.2. SB anomaly 

As described in SI our image subtractions degrade 
with increasing galaxy surface brightness, leading to in¬ 
creased flux scatter (see Figsures 1^ and [T0)l . The ori¬ 
gin of this SB anomaly has not been identified, but we 
speculate that it may be caused by an underestimate of 
the pixel flux errors in resampled images in the vicinity 
of bright galaxies. In particular, resampling introduces 
pixel-to-pixel correlations in the galaxy profile which are 
not included in our estimate of the PSF-fitted uncer¬ 
tainties. Other possibilities include subtle problems in 
the astrometric solution, the PSF determination, or the 
coadding of exposures. 

To check for the possibility that we introduced the SB 
anomaly in our customized version of hotPants, we have 
run a few tests using the publicly available version. We 
find that the subtracted images look very similar, and 
that our version results in notably fewer outlier fluxes. 
We are therefore confident that we have not introduced 
bugs to cause the SB anomaly. 

In the literature on transient-search pipelines we could 
not find a quantitive analysis on the effect of subtrac¬ 
tions on bright galaxies. However, there are some inter¬ 
esting clues in the final-photometry results reported by 
Pan-STARRSl and SNLS. In the recent Pan-STARRSl 
cosmology analysis, which uses the same underlying sub¬ 
traction technique as our Diffimg, their light curve fits 
have a reduced distribut ion with a larger h igh-side tail 
than expected (see Fig 6 in iRest et al.l (j2014[ B. They at¬ 
tribute this effect to subtraction artifacts on bright galax¬ 
ies, which is similar to our SB a nomaly. _ 

In the SNLS final photometry (jAstier et al.ll20I^ , they 
use a scene modeling techni que with stacked imag es, orig¬ 
inally developed for SDSS (jHoltzman et al.ll2008 ), which 
does not use resampled images. As a function of total 
SN + galaxy brightness, they find no evidence f or flu x 
bias or scatter (see Figs. 7 and 10 in lAstier et al.l (I2013D . 
which is encouraging that the SB anomaly can be re¬ 
solved in the offline analysis. It is not clear if their lack 

See Table 2 in PIO for SNLS exposure time in each band. 


of SB anomaly is due to a different photometry method, 
their astrometric precision being an order of magnitude 
better compared to our Difflmg|3 or because they do 
not probe sufficiently bright galaxies to see the effect. 

We are actively developing a final-analysis photome ¬ 
try method similar to that in iHoltzman et al.l (j2008li : 
lAstier et al.l (|2013D . but the SB anomaly may not get 
resolved for finding transients with Diffimg. The SB 
anomaly’s impact on discovering SNe la, however, is 
quite limited because of their brightness at low redshifts 
where the SB anomaly is most pronounced, and because 
only 2 detections are needed among of the many above¬ 
threshold observations. The main impact is that the 
larger flux uncertainties at low redshift slightly degrade 
the classification performance of the PSNID program. 

In contrast to bright SNe la, the SB anomaly can 
have a more dramatic effect on detecting and measur¬ 
ing fluxes for faint or fast transients, such as CC SNe 
or kilonovae. For example, kilonova models for neutron- 
star (NS) mergers suggest optical signals that are much 
dimmer redder, and s hort-lived compared to SNe la 
(|Barnes fc KaserJ l2013f ). Using Diffimg to search for 
such events in very nearby galaxies, the SB anomaly 
could significantly degrade the detection efficiency, and 
those that are detected could have color uncertainties 
much larger than expected from photo-electron statis¬ 
tics, thereby making it difficult to distinguish kilonovae 
from other astrophysical transients. 

To further diagnose the SB anomaly. Fig. [14] shows the 
autoScan score distribution for i band fake detections in 
the two deep fields (X3,C3). AutoScan assigns a score 
near zero to a clear artifact, and a score near one to a 
cleanly subtracted point-source transient; scores above 
0.5 are used to make candidates. The upper-left panel 
in Fig. [T4| shows the autoScan score distribution for all 
of the i band detections; this reference distribution is 
strongly peaked near one, showing that most of the de¬ 
tections are from good subtractions. The remaining pan¬ 
els show the autoScan score distribution, in bins of tosb, 
for the small subset of > 3(T flux outliers. For the bright¬ 
est SB range (20 < msB < 21) the autoScan scores are 
all close to zero, indicating that these are visibly poor 
subtractions. As the SB decreases, the autoScan scores 
improve. We have checked the distributions of PSF, sky 
noise and ZP, and find no significant difference between 
the outliers and the reference; hence there is no apparent 
correlation of the SB anomaly with observing conditions. 

For PSNID light curve fitting we could remove the few 
observations that fail autoScan but we do not currently 
have the infrastructure to apply this requirement to the 
many non-detections that are often more numerous than 
the autoScan failures. As described in [JS] we have chosen 
instead to model the increased flux scatter and inflate the 
flux uncertainties based on msB- 

Finally, we note that our characterization of the bright- 
galaxy subtraction artifact is a dependence on a single 
parameter: msB- While this description is adequate 
to classify newly discovered SNe for spectroscopic ob¬ 
servations, a more accurate description may be needed 
for dimmer transients (e.g., kilonovae), and the Hubble- 

While the SNLS final-photometry pipeline has much better 
astrometric precision than our Diffimg, the SNLS search pipeline 
(PIO) and Diffimg have similar astrometric precision. 
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diagram analysis if this effect persists in the final pho¬ 
tometry. For example, the SB anomaly could also depend 
on the exposure conditions and the SB gradient at the 
SN location. 
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Fig. 14.— autoScan score distribution for fakes in the two deep 
fields, i band. Score is 0 for a deary bad subtraction and 1 for a 
good subtraction. Upper left panel shows the reference distribution 
for all fakes. Remaining panels probe the SB anomaly for a subset 
of fakes that are faint > 23), are detected by SExtractor, 

satisfy the cuts in Table 0 and have a measured flux more than 
3 (t away from the true flux. Each panel indicates a different msB 
range. 


7.3. Host-galaxy Matching 

The SN science analyses will rely mainly on photomet¬ 
ric classihcation, and the redshifts will com e from host 
galax y spectroscopy, primarily from OzDES (jYuan et al.l 
I 2 OI 5 II . The spectroscopic redshifts are very accurate in 
principle, if the correct galaxy is matched to each SN. We 
have used fakes to measure the SN-host matching perfor¬ 
mance in Yl, and found a 99% success rate. However, 
our fakes are preferentially distributed close to the galaxy 
cores with too few events in the disk tails, and thus the 
SN-host matching result from fakes is too optimistic. 

We are therefore preparing to test SN-host matching 
with an independent set of fake locations based on more 
realistic galaxy profiles from semi-analytic models that 
are fit to Sersic profiles. As we obtain more accurate 
DES galaxy profiles in future analyses, we will be able 
to use our own data to evaluate the SN-host matching 
efficiency. Also note that fake locations can be rapidly 
generated and analyzed at the catalog level since there is 
no need to overlay SN fluxes on images to process with 
DiffImg. The eventual goal is to update the simulation 
to include a model of outlier redshifts from mis-matched 
host galaxies. 


8. CONCLUSIONS 

We have assembled a pipeline capable of using hun¬ 
dreds of CPU cores to process up to 170 GB of raw 
imaging data in less than a day, with the goal of dis¬ 
covering astrophysical transients. Eor the subtracted im¬ 
ages produced by Diff Img in Yl, the typical number of 


SExtractor-detected objects per band is a few hundred 
thousand per 3 deg^ field, and the vast majority (> 90%) 
are subtraction artifacts. Selection requirements and au¬ 
tomated scanning reduce the artifact fraction down to 
25%, and ~ 10^ detections per band (Table 01). The 
number of detections per single-epoch visit is ^ 130 per 
deg^. 

The number of science candidates, requiring a detec¬ 
tion on 2 separate nights, is 1040 per deep field, and 680 
per shallow field (Table E]). Our MC simulation predicts 
that roughly 27% of the discovered transients are SNe la 
or CC SNe. Another ^ 30% are artifacts, and most of 
the remaining candidates are AGN or variable stars. 

We have implemented extensive monitoring in Diff Img 
based on overlaying fake SNe la near galaxies on the 
search images. Gomparing the Diff Img efficiency for 
fakes to the efficiency from MG simulations shows that 
the Diff Img performance is close to what is expected. 
The main defect of Diff Img is the SB anomaly in which 
larger host-galaxy surface brightness results in larger 
flux-scatter that is not described by the uncertainty (see 
Figures 0] and Uni). There are other small fake-MG dis¬ 
crepancies in the efficiency (e.g., Fig. 03; it is not clear 
if the cause is a more subtle Diff Img defect, or if the MG 
simulation is too optimistic. 

As a rigorous demonstration of our monitoring tech¬ 
nique, we performed a very preliminary photometric clas¬ 
sification analysis on real (non-fake) data, and compared 
the resulting Zphot distribution to a MG simulation. In¬ 
puts to the MG simulation include observed conditions 
(PSF, ZP, sky noise) and the Diff Img behavior measured 
with fakes (efficiency vs. S/N and anomalous flux scatter 
vs. SB). The resulting data-MG agreement is reasonable 
in both the deep and shallow fields (Fig. [13]) • 

Finally, the results presented here are based on fully re¬ 
processed data after the first two DES seasons. Diff Img 
issues during Yl and Y2 resulted in some poor subtrac¬ 
tions, but with recent Diff Img improvements and a reli¬ 
able model of the flux uncertainties, we expect our spec¬ 
troscopic target selection to be more efficient and more 
automated in the remaining seasons. 
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APPENDIX 

A. MC SIMULATION TO PREDICT THE Diffimg EFFICIENCY 

The fast MC simulation of SNe la is from SNANA (|Kessler et al.l[200^ . It uses the exact same generation parameters 
as those used to generate the fakes ( 113.2.21) : parent populations of color and stretch, intrinsic scatter model, and a 
random galaxy location in proportion to its sur face brightness density. For studies requiring SN CC we use the SNANA 
simulatio n as described in (iKessler et al.ll20inaD . For studies requiring t he absolute rate, we use the SN la volumetric 
rate from iDildav et aO (j2008h and the CC rate from iBazin et al.l (j2009f l. Each simulated epoch corresponds to a real 
observation in the survey where the model magnitude is converted to an equivalent forced-photometry flux using the 
measured ZP. The observed PSF and sky noise at each epoch are used to predict the measurement uncertainty, 

4iM = [F+iA-b-RMSl)] (AI) 

where ctsim is the uncertainty in photoelectrons, F is the flux, A = [27r J PSF^(r, 6)rdr]~^ is the noise-equivalent area, 
and b is the effective sky level including dark current, readout noise, and noise from the host galaxy. RMS a is an 
empirical error scaling of the sky noise that increases with the local surface brightness as shown in Figures [9] and HOl this 
term accounts for the SB anomaly: systematic subtraction problems near bright galaxies. While the measured RMSa 
curves are used to compute anomalous fluctuations in the measured fluxes, the reported uncertainties are computed 
with RMSa = I in the same way as the data. 

The simulation includes the candidate selection requirement of a detection on two separate nights. The detection 
efficiency is computed from the cs/n curves in Fig. [8l A simulated detection requires cs/n > r, where 0 < r < I is a 
random number. 

Finally, the simulated light curves are stored in data files and analyzed in exactly the same way as transients (fakes 
or real events) found by Diffimg. 


B. PHOTOMETRIC ANALYSIS AND SELECTION REQUIREMENTS 

Here we describe a photometric analysis and selection requirements to obtain a high-purity sample of SNe la in the 
first season of the DES-SN program. The goal of this analysis is to compare the Zphot distribution for data and the 
MC simulation. Using the SALT-II model, light curve fitting is done with the SNANA program snlc_fit.exe. For 
each candidate, the 5 fitted parameters are (I) time of peak brightness (to), (2) SALT-II color parameter (c), (3) 
SALT-II stretch parameter (xi), (4) SALT-II amplitude (xq), and (5) photometric redshift (^phot)- 
The first fit iteration chi-squared (y^) is computed in the usual manner: from the data-model flux-difference for 
each epoch, and the quadrature sum of the data and model uncertainties. Since the model uncertainty depends on the 
fitted parameter Zphot, the second fit iteration chi-squared (xf) is 

xl = xi + xl where x^ = XI ^ ■ (BI) 

e 

The index e is the epoch index and trf is the quadrature sum of the data and model-uncertainty from the first fit 
iteration in which there is no Xa term. While the af add an irrelevant constant to xi, it has the effect of making Xa 
small. The analysis selection requirements are as follows: 

1. three bands with at least one observation satisfying S/N> 5. 

2. at least I observation with T^est < —2 days, where T^est = 7obs/(l + -^phot)- 
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3. at least 1 observation with T^est > +10 days. 

4. SALT-II stretch parameter \xi\ <4 

5. 0.02 < Zphot < 2 

6 . fit probability Pfn > 0.1, calculated from fit y^/dof. 

7. Ix^l < 2.5. 

8 . NN requirement described below. 


The NN analysis is based on the four-dimensional space of xi, c, 2phot and rriB- The first three variables are from 
the SALT-II light curve fit (see above), ms is the true rest-frame i?-band magnitude as described in Sec 4.3 of 
iKessler et al.1 (j2013ll . and is not the naive best-fit model magnitude. For a given set of fitted parameters, the NNs are 
simulated events that satisfy a four-dimensional distance constraint. 




(c-c')^ (xi-x()2 (^phot - 2phot)^ 


(B2) 


where the primed quantities are the fitted parameters from a simulated training sample that includes SNe la and CC 
SNe events. The optimal distance-metric parameters (A^, A^,,^, A^, A^) are trained with the simulation to maximize 
the product of the SN la purity and the efficiency. The final selection requirement is that for simulated neighbors 
satisfying Eq. IB21 more than half are true SNe la with at least Icr confidence. 
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